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Abstract 

We investigate the sensitivity of several observables to the density dependence of the symmetry 



potential within the microscopic transport model UrQMD (ultrarelativistic quantum molecular 



^ . dynamics model). The same systems are used to probe the symmetry potential at both low and 

h ; 

high densities. The influence of the symmetry potentials on the yields of tt , tt + , the it /tt + ratio, 
the n/p ratio of free nucleons and the £/ 3 He ratio are studied for neutron-rich heavy ion collisions 
( 208 Pb+ 208 Pb, 132 Sn+ 12i Sn, 9e Zr+ 9G Zr) at E b = 0.4AGeV. We find that these multiple probes 
provides comprehensive information on the density dependence of the symmetry potential. 

PACS numbers: 24.10.Lx, 25.75.Dw, 25.75.-q 



* Fellow of the Alexander von Humboldt Foundation. 



1 



The nuclear equation of state (EoS) for isospin-asymmetric nuclear systems has recently 



where 6 = (p n — p p )/(p n + p p ) is the isospin asymmetry defined by the neutron (p n ) and 
proton (p p ) densities. Cq is the energy per nucleon for symmetric nuclear matter and e sym (p) 
is the bulk symmetry energy. The nuclear symmetry energy term e sym (p) is very important 
for the understanding of many interesting astrophysical phenomena, but it is also plagued 
by large uncertainties. An extreme variation of the EoS is observed for neutron-rich matter 
by using different versions of the Skyrme interaction as well as other non-relativistic effective 
interactions |2j. The EoS of neutron-rich matter can be studied using a density-dependent 
relativistic hadron mean-field theory like the DDH3p5 and DDHp models P|. The density 
dependence of the symmetry potential deduced from relativistic and non-relativistic interac- 
tions might be different. Therefore, a comparison of the isospin effect in isospin-asymmetric 
heavy ion collisions (HICs) with both relativistic and non-relativistic effective interactions is 
important. Furthermore, most of the studies on sensitive probes of the density dependence 
of the symmetry potential are concerned with the behavior in the low-density or in the high- 
density region separately. Thus the simultaneous study of observables, which are sensitive 
to both the low-density and the high-density part of symmetry potential, in one reaction 
system is necessary. This allows to extract the full information on the density dependence 
of the symmetry potential which relates to specific effective interactions. 

To study the isospin effects in intermediate energy HICs, heavy stable nuclei with an 
isospin asymmetry as large as possible should be used. They have the advantage that both, 
the availability and the beam intensities are much higher than for radioactive beams and, 
hence, high statistics experimental data can be obtained. This is necessary to study the 
relatively subtle effects of different symmetry potentials. The reaction 208 Pb + 208 Pb, for 
example, has an isospin-asymmetry 5 ~ 0.212, i.e., a neutron to proton ratio N/Z ~ 1.54, 
which is already quite large 0]. 

In a previous work, we have adopted the following forms of the symmetry potential 




e(p, 5) = e (p, 0) + e sym (p)<5 2 



(1) 
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to mimic the strong variation of the density dependence of the symmetry potential at high 
densities as given by the Skyrme-type interactions. Here, u = p/po is the reduced density 
and a is the reduced critical density. When u > a, the symmetry potential energy is negative 
in F 2 . We choose 7 = 1.5 (a = 3) for i*\ (F 2 ) and name it F15 (Fa3). In the present work we 
further adopt two other forms of the symmetry potential based on the relativistic mean field 
theory, namely DDH3p<5* and DDHp*. Both are based on the extended QHD model (see 
Ref. 3j). The density dependence of the symmetry potentials for the DDH3p5 and DDHp 
models is due to the density dependence of the coupling constants (g CT jv, gwN, S p n, and g^). 
Fig. H shows the density dependences of these symmetry potentials as well as the linear case. 
For comparison, we fix the symmetry potential energy at normal nuclear density for all four 
forms of the symmetry potential. The symmetry energy coefficient is adopted to be So = 34 
MeV js, E|. It is seen from Fig. QJthat, for u < 1, DDH3p<5* is very close to Fa3, both of 
them lie between DDHp* and F15. For u > 1, they show a rather different behavior and 
approach F15 and DDHp*, respectively. 

In the present work, central collisions are studied for the three neutron-rich systems 
208 Pb + 208 Pb (5 ~ 0.212), 132 Sn + 124 Sn {5 ~ 0.218), and 9e Zr + 96 Zr {5 ~ 0.167) at 
Efy = OAA GeV. The ultrarelativistic quantum molecular dynamics model (UrQMD, version 
1-3) jsllil 10, 11 1 is adopted for the calculations including the Coulomb potential for charged 
pions. A 'hard' Skyrme-type EoS (K = 300 MeV) without momentum dependence is used 
for the calculations. It should be noted that the uncertainty of the isospin independent EoS 
is still large (see Refs. jfllisj]), and the contribution of the momentum dependence in EoS 
is also important to the dynamics of HICs (see Ref. Q]), both of them influence the results 
but do not change the effect of the density dependence of symmetry potential on the ratios 
of different charged particles strongly [l^ . 

Figs. El and H3 show the rapidity and the transverse momentum distributions of 7r~ and 
7T+ for central (b = - 2 fm) 208 P6 + 208 Pb collisions at OAA GeV as calculated with the 
symmetry potentials F15 and Fa3. Fig. El also shows the corresponding 7r~ /tt + ratios. 
While the ir~ yields are quite sensitive to the symmetry potential at mid-rapidity and at 
low p£ m , the 7r + do not exhibit this sensitivity due to the neutron-rich environment jl^ . 
The calculations for the 132 Sn + 124 Sn and 96 Zr + 96 Zr reactions are similar to those in 
Figs. 121 and El To analyze the differences in the 7r~ and tt + yields and the 7r _ /tt + ratios (as 
calculated with the four different forms of the symmetry potential) more quantitatively, we 
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FIG. 1: Parametrizations of the nuclear symmetry potential DDHp*, DDH3p<5*, F15, Fa3, and the 
linear one as a function of the reduced density u. 

define a variable = (JQ — Xj)/Xj, where the subscripts "z" and "j" denote different 
forms of the symmetry potential. "X" represents the calculated quantity, such as the pion 
yield or the tt~/tt + ratio. Here, % denotes the symmetry potential that yields a larger X 
value than the one with the symmetry potential j, i.e., we always choose Xi > Xj. Table U 
(for 208 Pb + 208 Pb, titled as "P6208") and Table El ( for 132 Sn + 124 Sn, titled as "Snl32") 
give the results of Dij for different conditions and physical cuts: the subscript "tot" indicates 
the results without any kinematical cut, "y ~ 0" denotes particles emitted at mid-rapidity 
(y c in the center-of-mass system), and the subscript "p t ~ 0.1" means the particles with 
transverse momenta p£ m ~ O.lAGeV /c . Here, the index 1 represents the "DDH3p5*" 
symmetry potential, 2 represents "DDHp*", 3 represents "F15", and 4 is the "Fa3" case. 
Dmax provides the maximum value of all Dy. The differences are maximal for the D 2 3 
values, these are the differences between the F15 and DDHp* symmetry potentials, which 
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also exhibit the largest deviations in the density dependence of the symmetry potential at 
densities u ~ 1 — 2.6 (see Fig. QJ). Tables HI and ITT1 demonstrate that the values of Dij 
for 7r + are always smaller than those for 7r~, this is consistent with Figs. 121 and El By 
comparison, we find that the ir~ /ir + ratio at p£ m ~ 0.1 GeV/c is most suitable to probe the 
symmetry potential because of the larger ir~ and 7r + yields in this transverse momentum 
region. We divide all into two groups in Tables U and |Hj In general, the values in the 
first group of D 2 i, -D 43 and D max {= D23) are much larger than those in the second group 
of -D13, -D41, D24. For the former case the values are generally larger or close to 10 except 
for the 7r + , while for the latter one they are all smaller than 10, except D13 for the n~ /n + 
ratio at p£ m ~ O.lGeV/c with a value of 10. Surprisingly, _D 41 lies in the second group, 
although Fa3 corresponds to a soft symmetry potential and DDH3p5* corresponds to a stiff 
symmetry potential. The reason can be found in Fig. ^ the difference between these two 
symmetry potentials is not large for reduced densities u = 1 ~ 2 where most of the pions 
are produced Q|. It means that more accurate measurements are needed in order to probe 
the density dependence of the symmetry potential with pions, except for the extreme cases. 
Furthermore, from Tables HI and HTl we observe Xddh p * > Xpad, > Xddh3 P 8* > Xpi5- This 
order agrees with the order of the four forms of the symmetry potentials at baryon densities 
u — 1 — 2 (see Fig. Baryon densities of u ~ 1 — 2 are the most relevant ones for pion 



production as we have shown in 



151 ] . Comparing the values shown in Tables U and |Hj 



we see that the values are a little smaller for 208 p5-|_ 208 p5 as compared to 132 Sn+ 12A Sn 
with radioactive beams. 

Probes being sensitive to the density dependence of the symmetry potential at subnormal 
densities have been studied extensively in recent years. Usually a form (p/po) 7 is assumed 
for the density dependence of the symmetry potential ^(|. From Fig. ^ it follows that the 
behavior of the symmetry potential at subnormal and high densities are quite different. To 
gain information on the symmetry potential in both, the high and the subnormal density 
regions simultaneously, observables from the same reaction system are needed, which are 
sensitive to both density regions. Most of the free nucleons, intermediate mass fragments 
and light charged particles are emitted in the expansion stage, that is, at subnormal densities. 
Therefore, we further study the ratios of free neutrons and protons, the isotope distributions 
of light and intermediate mass fragments, as well as the ratio of emitted tritons and 3 He 
which is also proposed as a sensitive probe for the symmetry potential at subnormal densities 
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FIG. 2: The rapidity distributions of ir and tt + for central 208 Pb + 208 Pb collisions with different 
symmetry potentials (see text). 

\v\ . A conventional phase-space coalescence model is used to construct the clusters, 
in which the nucleons with relative momenta smaller than P and relative distances smaller 
than i?o are considered to belong to one cluster. In this work Pq and Rq are chosen to be 
0.3GeV/c and 3.5 fm, respectively. Another set, Pq = 0.25GeV/c and Ro = 3.0 fm, alters 
the yields of free nucleons and clusters strongly, but it does not change the dependence of 
the t/ 3 He ratios on the symmetry potentials. The freeze-out time is taken to be 150 fm/c. 

n n 

In Ref. pj| a modified coalescence model [19] was adopted. It was found that a modification 
of the coalescence criteria changes the yields of light charged particles but it does not change 
the i/ 3 He ratio. 

Fig. H] shows the transverse momentum distribution of the n/p ratios of emitted unbound 
nucleons as calculated with the four different forms of the symmetry potential. The n/p ratio 
obviously depends strongly on the choice of the symmetry potential. In the low transverse 
momentum region the n/p ratio is the largest for DDHp* and the smallest for F15. Obviously, 
nucleons with low transverse momenta are mainly emitted from the low density region. 
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p™ [GeV/c] 

FIG. 3: Transverse momentum distributions of 7r - and 7r + from central 208 p5_j- 208 p5 collisions at 
Eh = OAA GeV for different symmetry potentials. The tt~/it + ratios are also shown as a function 
of transverse momentum. 

TABLE I: Values Dij of the total 7r ± yields (at mid-rapidity or p£ m ~ O.lGeV/c) and of the 
corresponding ir~ /ir + ratios for various combinations of the density-dependent symmetry potentials 
for 208 Pb + 208 Pb reactions (see text). 
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TABLE II: The same as in Table [I] but for 132 Sn + 124 Sn reactions. 
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From Fig. ^ it follows that the difference between Fa3 and DDH3p5* is small at densities 
p < po- Correspondingly, the n/p ratios calculated with Fa3 and DDH3p<5* are also very 
close at lower transverse momenta. For emitted nucleons with transverse momenta larger 
than ~ 700 MeV/c, the n/p ratios in the Fa3 and DDH3p<5* cases are close to those in the 
DDHp* and F15 cases, respectively. This reflects the fact that E£°^ for Fa3 (DDH3p<5*) 
is close to DDHp* (F15) when p > p . Free nucleons with transverse momenta larger 
than ~ 700MeV/c are mainly squeezed out from higher densities j^. Thus, the different 
behavior of the n/p ratio as a function of transverse momentum goes in line with the behavior 
of the symmetry potential at low and high densities as shown in Fig.d The relative difference 
between the n/p ratios for transverse momenta larger than ~ 700 MeV/c as calculated with 
the soft and stiff symmetry potentials is on the same order as for the it~/tt + ratios. 

To see the different emission times for nucleons with high and low transverse momenta, 
Fig. El shows the transverse momentum distributions of free neutrons and protons (upper 
plot) as well as the corresponding n/p ratios (lower plot), emitted from the mid-rapidity 
region (|(y^| = \y c /ybeam\ < 0.3) at reaction times t = 50 and 150 fm/c. In the upper 
plot of Fig. El only the results for F15 are shown while the ratios for both the F15 and 
DDHp* symmetry potentials are compared in the lower plot. Nucleons with large transverse 
momenta are emitted in the early stage of the reaction and nucleons with small transverse 
momenta are emitted at later time. Furthermore, we see an obvious change of the transverse 
momentum distribution of the n/p ratios with time, especially in the low transverse momen- 
tum region, in particular in the case of the soft symmetry potential DDHp*. This can easily 
be understood from Fig. ^ in the cases of the soft symmetry potentials DDHp* and Fa3, 
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FIG. 4: Transverse momentum distributions of free neutron/proton ratios for different density- 
dependent symmetry potentials. 

neutrons are strongly attracted into the high density region by the symmetry potential at 
the compression stage. Within the expansion period, neutrons are driven to lower densities, 
eventually leading to the subnormal densities. At subnormal densities the soft symmetry 
potential energy is higher than the stiff one. Hence, more neutrons with lower transverse 
momenta are emitted for the cases with soft symmetry potentials DDHp* and Fa3 at the 
later stage (at 150 fm/c) as compared to the cases with the stiff symmetry potentials. For 
the stiff symmetry potential case, there is no significant change of the transverse momen- 
tum distributions of the n/p ratio from t = 50 fm/c to 150 fm/c. In the high transverse 
momentum region, the n/p ratios for both cases decrease with time during the expansion 
phase, which happens mainly due to the extra Coulomb effect on the protons. 

However, experimentally it is difficult to measure the spectra of neutrons with high pre- 
cision. Therefore, the production of light charged particles with large but different N/Z 
ratios as well as isotope distributions for some intermediate mass fragments could be more 



useful Oil 



231 • Fig. El shows the isotope distributions of 1_3 H, 3 " 6 He, 6_8 Li, 8 - n Be, 10 - 13 B 
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FIG. 5: Transverse momentum distributions of free neutrons and protons (upper plot) as well 
as the corresponding n/p ratios (lower plot) at mid-rapidity (|(yc°^| < 0.3) and at reaction times 
t = 50 and 150 fm/c. The F15 and DDHp* symmetry potentials are chosen. The arrows in the 
upper plot show the maxima of the distributions. 
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FIG. 6: Isotope distributions of H, He, Li, Be, B and C with the symmetry potentials F15 and Fa3 
(see text). 

and 12 ~ 15 C in central Pb+Pb collisions as calculated with the symmetry potentials F15 and 
Fa3, respectively. In general, the differences between the isotope yields for the F15 and Fa3 
cases are smaller than the differences in the n/p ratio of emitted nucleons. However, as the 
isotopes become more neutron- rich, the corresponding yields decrease fast. Secondly, with 
an increase of the neutron number in the light isotopes of H, He and Li, the differences 
between the yields calculated with F15 and Fa3 enlarge (diminish) for symmetric (asym- 
metric) isotopes, while for the heavier Be, B, and C elements the isotope distributions are 
almost the same for F15 and Fa3. This can be understood based on the density dependence 
of the symmetry potentials shown in Fig. ^ for He and Li isotopes, a weak enhancement 
of isospin symmetric isotopes with Fa3 is shown but this is not observed for neutron rich 
isotopes. This is because the n/p ratio for free nucleons, as calculated with Fa3, is much 
larger than that with F15, which leads to the reduction of the N/Z ratios for H, He, and Li 
when calculated with Fa3. 

The light charged particles like triton and 3 He ought to be particularly noted, because they 
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exhibit a considerable difference in the isospin asymmetry, as well as large yiel ds in heavy ion 
collisions at SIS energy, E b < 1A GeV. Both the FOPI IgJ and the EOS [gj collaborations 
have found that - in central Au+Au collisions - 3 He clusters have a systematically larger 
mean kinetic energy than tritons up to bombarding energies of 400^4 MeV. This trend is 
also found in our calculations. Fig. |U] shows the yield of tritons is larger than that of 
3 He. To see the symmetry potential effect more clearly, Fig. [7| shows the t/ 3 He ratios as 
calculated with the four different forms of the symmetry potential used in this work. The 
dashed line is for the case with the rapidity cut and the solid line is for the case without 
any kinematical cut. Fig. shows that the t/ 3 He ratios are more strongly enhanced with 
a softer symmetry potential at mid-rapidity \y^'\ < 0.5 than without any cut. This is 
because neutron-richer matter is formed in the mid-rapidity region with an increasingly softer 
symmetry potential. Thus, more (less) tritons ( 3 He) are produced with a softer symmetry 
potential as compared to the case without any cut. The magnitudes of the calculated t/ 3 He 
ratios can be ordered as follows: Xbdh p * > Xddh3 P s* > Xp a3 > Xpi§. This is a little 
different from the ordering of the 7r~ /h + ratios. However, the difference of t/ 3 He is very 
small between Fa3 and DDH3p5*. D max is about 22% for the t/ 3 He ratio in full phase space 
(47r) and D max ~ 27% for the \(y^\ < 0.5 case. This is larger than in the case of pions 
as given in Tables |l] and |Hj Comparing with Fig. ^ it becomes clear that the ordering of 
the t/ 3 He ratios for the four different forms of the symmetry potential corresponds to the 
ordering of the E^°^ at subnormal densities! This holds, regardless of the behavior at u > 1. 
Therefore, we consider the t/ 3 He ratio as a sensitive probe of the density dependence of the 
symmetry potential at subnormal densities. 

In summary, we have investigated the influence of different density dependences of the 
symmetry potential on charged pion yields, the 7t~/tt + ratio, the n/p ratio for free nucle- 
ons, the isotope distribution of H, He, Li, Be, B, C and the t/ 3 He ratio for 208 Pb + 208 Pb 
(5 ~ 0.212), 132 Sn + 124 Sn (5 ~ 0.218), and 96 Zr + 96 Zr (5 ~ 0.167) at E b = 0.4AGeV, 
respectively. We have shown that the negatively charged pion yields, and the it~ /tt + ratios 
at low momenta are very sensitive to the density dependence of the symmetry potential at 
high densities, while the n/p ratio, i.e., for free nucleons, and the t/ 3 He ratio are sensitive 
to the potential at subnormal densities. The n/p ratios at high transverse momenta show 
sensitivity to the behavior of the symmetry potential at high densities. To obtain the full 
information on the density dependence of the symmetry potential, several sensitive probes 
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FIG. 7: The ratios of the triton and 3 He yields with and without rapidity cut as calculated with 
four different forms of symmetry potentials. 

at both high and low densities were analyzed simultaneously. The comparison of the results 
of different systems shows that it is suitable to use stable but highly isospin-asymmetric 
systems to study the effects of different symmetry potentials. 

Here, it should be mentioned again that concrete values for the yields of emitted pions, 
neutrons, protons, tritons, and 3 He may be influenced by the uncertainty of the isospin- 
independent (and the momentum dependent) part of the EoS. However, we expect that the 
dependence of the n/p , tt~ /tt + , and t/ 3 He ratios on the symmetry potential is not strongly 
affected. 
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